Critical velocity for the vortex core reversal in perpendicular bias magnetic field 
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For a circular magnetic nanodot in a vortex ground state we study how the critical velocity v c of 
the vortex core reversal depends on the magnitude H of a bias magnetic field applied perpendicularly 
to the dot plane. We find that, similarly to the case H = 0, the critical velocity does not depend 
on the size of the dot. The critical velocity is dramatically reduced when the negative (i.e. opposite 
to the vortex core direction) bias field approaches the value, at which a static core reversal takes 
place. A simple analytical model shows good agreement with our numerical result. 

PACS numbers: 75.78.Fg, 75.60.Jk, 75.78.Cd, 75.75.-c 
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A magnetic vortex is a curling magnetization distribu- 
tion in flat magnetic submicron dots, with the magneti- 
zation pointing perpendicularly to the dot plane within 
the ten nanometer size vortex core. The vortex ground 
state corresponds to a deep energy minimum when the 
dot lateral sizes fit the conditions for vortex stability [l[ . 
This unique magnetic object has attracted much atten- 
tion recently because of the fundamental interest to spe- 
cific properties of such a nanoscale spin structure. The di- 
rection of the core polarization ('up' or 'down') can store 
a bit of information, and this is of considerable practical 
interest for applications in magnetic memory technology. 
Several approaches can be used to switch the bit, i.e., 
to reverse the vortex core. It has been shown that a 
static magnetic field can reverse the core if its magni- 
tude reaches sufficiently large values, typically of several 
kOe [2, l3j ■ The reason for that is a large energy barrier 
between the vortex states with 'up' or 'down' core polar- 
izations. Alternately, the magnetic core can be switched 
at zero static magnetic field, if it is excited by a vari- 
able (oscillating or pulsed) in-plane field or by a spin- 
polarized current 0, H, S 0] ■ The reversal occurs if the 
core velocity reaches a certain critical value v c , which is 
defined solely by the magnetic parameters Q. Very re- 
cently, the vortex core switching has been observed at 
intermediate experimental conditions in Ref. 0. In this 
work, a static perpendicular magnetic field together with 
a small oscillating in-plane field was applied to a nan- 
odot in a vortex state. The frequency of the excitation 
was swept and the resonant vortex motion was detected. 
At a given magnitude of the exciting field, the vortex core 
was switched when the static bias magnetic field reached 
a critical value. 

In our numerical study, we calculate the critical ve- 
locity of vortex core reversal as a function of a static 
out-of-the plane magnetic field. We find that this criti- 
cal velocity, similarly to the case of zero applied field, is 



independent of the dot sizes, but depends on the mag- 
netic parameters of the dot material. The critical velocity 
drops significantly with the increase of the magnitude of 
the negative (opposite to the vortex core direction) bias 
magnetic field. 

We consider a circular magnetic nanodot in a vortex 
ground state. A static bias magnetic field is applied per- 
pendicularly to the dot plane (along the z-axis). It is 
considered to be positive if parallel (and negative if an- 
tiparallel) to the initial direction of the vortex core (or 
vortex core polarization), see the inset to Fig. [TJ The 
vortex motion is excited by a d.c. spin polarized cur- 
rent flowing perpendicular to the dot plane 10]. The 
spin polarization of the current is along the z-axis. The 
vortex motion is calculated numerically using our micro- 
magnetic code Two dots with a diameter 2R = 300 
nm and thicknesses w = 20 nm and 30 nm were con- 
sidered. Magnetic parameters mimic those for NiMnSb 
used in Ref. [9j : the saturation magnetization is 4ttM s = 
0.69 T, the easy plane anisotropy field H a = 0.185 T, the 
exchange stiffness A = 1 x 10~ n J/m and the Gilbert 
damping a = 0.01. The mesh cell size is 1.5 xl.5 x 5 
nm 3 . 

Due to the excitation by the spin current, the vortex 
core starts to gyrate with gradually increasing radius. 
Correspondingly, the core velocity is increasing until it 
eventually switches. Prior to the core reversal, a region 
with negative values of M z component (a 'magnetic dip') 
is formed at the inner part of the core trajectory. When 
the core velocity reaches the critical value, this dip splits 
into a vortex with negative polarization and an antivor- 
tex [6] . The antivortex annihilates with the original vor- 
tex core, and in the end only the vortex with a negative 
polarization remains. After the core reversal the spin- 
polarized current starts to damp the core gyration, thus 
slowing the core motion and, eventually, bringing the re- 
versed core to the equilibrium position in the dot center 
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The critical velocity v c is determined as the maximum 
core velocity on the gyration trajectory. This maximum 
is reached just before the core switching. The core ve- 
locity is calculated as the time derivative of the core po- 
sition, which in its turn is extracted from the magneti- 
zation distributions printed out each 0.05 ns. At fields 
larger than 0.12 T for the dot with w = 30 nm (corre- 
spondingly, larger than —0.04 T for the dot with w — 
20 nm) the core is expelled from the dot prior to the 
switching. At fields smaller than —0.55 T the mesh we 
use becomes insufficiently fine to be able to calculate the 
vortex dynamics accurately. The static switching field 
H r is equal to —5.9 kOe for the both dot thicknesses 
[121 ]. This approximate numerical value for H c was used 
below in our analytical calculations of the dependence of 
the critical velocity on the perpendicular bias field. The 
results of this calculation are shown by a solid line in Fig. 
1. 

The simulation results for the two dots are summarized 
in Fig. [T] (symbols). The critical velocity at H = is v c 
— 360 m/s for both dots. v c (H) increases for increasing 
positive fields (v c = 460 m/s for w = 30 nm at H = 0.12 
T). However it diminishes significantly for negative fields 
(v c = 40 m/s at H = - 0.55 T for both the dots). At 
moderate fields (\H\ < 2 kOe) v c scales linearly with H, 
and the slope is approximately 670 m(sT)~ 1 . 
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FIG. 1: (color online) Symbols: critical velocity v c as a func- 
tion of the magnitude of the magnetic field applied perpen- 
dicularly to the dot plane for dots with thickness w = 20 nm 
and w — 30 nm. Solid line: analytical prediction by Eq. [T] 

We find that for all the field values, when they can be 
calculated, the critical velocities determined for the two 
dots coincide. This fact is rather nontrivial as, owing 
to different thicknesses, many parameters of the vortex 
(such as the profile of t he p otential energy VF(X), where 
X is the core position [13j ; the separation of the vortex 
from the dot center and edges at the switching; the shape 



of the core) are very different for the two dots and depend 
differently on the field. From these results we conclude 
that, similarly to the case H = 0, for non-zero H, v c 
depends o nly on local properties of the vortex core spin 
structure [14|. 

The zero-field value of the critical velocity for both 
dot thicknesses, v c = 360 m/s, is in perfect agree- 
ment with the analytical prediction of the works of 
Ref. [H, [l5| taking into account the easy-plane 
anisotropy constant K — M s Ha/2 of NiMnSb: v c (0) = 
1.66M sA /27rA/(27rM s 2 + K), which gives v c = 340 m/s. 
In the following, we investigate the underlying physics 
responsible for the v c {H) behavior presented in Fig. 1. 
The vortex core dynamic reversal, as it was shown in Ref. 
[8j], originates from the self- induced dynamic gyrotropic 
field or gyrofield. This field is induced by the vortex mo- 
tion and its amplitude is proportional to the ratio v/p, 
where v is the velocity of the moving vortex and p is the 
vortex core radius. When the gyrofield reaches a critical 
value Hg T oc v c / p, the vortex core very rapidly reverses. 

We study how p scales with H for the two dot thick- 
nesses. We analyze magnetization distribution profiles 
for a static vortex in equilibrium at different fields to ex- 
tract the dependence p(H) [l^ . We find that although 
p(H) is different for the two dots (e.g., at H = 0, p = 20 
nm for the dot with w = 30 nm and correspondingly p 
= 18 nm for w — 20 nm), to the precision of our calcula- 
tion, p(H)/p(0) coincides for the two thicknesses, as can 
be seen on Fig. [2] From this result we conclude that the 
critical value of the gyrofield Hg T oc v c /p scales equally 
with the field H for the two dots. It also indicates that 
at non-zero external field the critical velocity v c (H) relies 
on the same vortex core reversal mechanism than at zero 
field, i.e. it is mainly determined by a competition of 
the gyrotropic and exchange fields within the core. The 
gyrofield deforms the core magnetization profile, whereas 
the exchange field tries to create a more uniform magne- 
tization distribution suppressing the core deformation. 

As can be seen from Fig. f2J the slopes of v c and p 
as functions of the perpendicular magnetic field are dif- 
ferent; indeed, v c (H) decreases noticeably more rapidly 
than p{H) at negative H . This means that the criti- 
cal value of the gyrofield H^ r decreases with negative H . 
This feature can be attributed to the fact that the de- 
formation of the core by the perpendicular bias field H 
leads to a decrease of the effective potential barrier that 
the gyrofield has to surpass to induce the vortex core re- 
versal. Therefore, the bias field provides two different 
actions on the vortex which help to switch the core. One 
is that, at a given core velocity v, the amplitude of the 
gyrofield H g increases with negative H as long as p is 
reduced. Second is that the critical value of the gyrofield 
Hg T that is required to switch the core becomes smaller 
at higher negative fields. For any field H, the vortex 
core reversal mode is an axially asymmetric mode like it 
was found for H = 0, Q . That is very different from 
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FIG. 2: (color online) Symbols: radius of the vortex core p 
as a function of magnetic field H, for two dots. Solid line: 
critical velocity for the dot with w = 30 nm. 



the axially symmetric reversal path involving the Bloch 
joint (BP) found in the simulations for static reversal 
3j . We also see this axially symmetric mode and the BP 
formation in our simulations of the static core reversal. 
But this axially symmetric BP mechanism is an ideal- 
ization, which leads to higher values of H c . It can be 
not realized practically due to unavoidable spontaneous 
symmetry breaking in real systems, e.g. induced by the 
thermal fluctuations. This is the vortex gyrotropic mode 
with a finite X that is responsible for the axial symmetry 
breaking. That is why the critical velocity v c (H) of the 
moving vortex is important. 

We can get a simple analytical expression for the 
v c (H). Let us consider a dot with static switching field 
H c . The value of the core radius at this field p c {H c ) is 
finite. The physical sense of p c (H c ) is the following: the 
vortex with positive polarization becomes unstable in the 
point H — H c when decreasing H. From the other side, 
it is reasonable to assume that the dependence v c (H) 
goes to when H approaches the static core reversal 
field H c ; i.e., we can assume that v c {H) is proportional 
to (1 — H/H c ) near H c . That immediately leads to the 
dependence: 



v c (H)=v c (0)(l-H/H c ), 



(1) 



The static field reversal and dynamic reversal mecha- 
nisms help each other leading to descending dependence 
of v c (H). Thus, the analytically estimated slope of the 
dependence v c {H) is dv c /dH — —v c (0)/H c = 610 m/s 
T (shown as a solid line in Fig. 1), that is close to the 
numerically simulated slope of 670 m/s T. These spec- 
ulations explain the main features of our simulations of 
v c (H), p c (H) presented in Fig. 1 and 2. 

In summary, our numerical study has demonstrated 
that there are two contributions to the process of the 
vortex core reversal in a magnetic dot subjected to a per- 
pendicular bias magnetic field: the static reversal mech- 



anism related to the instability of the vortex core with 
polarization directed against the bias field and the dy- 
namical reversal mechanism related to the vortex core 
deformation. While the first mechanism keeps the ax- 
ial symmetry of the vortex magnetization distribution, 
the second one breaks this axial symmetry and creates 
an "easy" core reversal path. Thus, the perpendicular 
bias magnetic field applied oppositely to the vortex core 
direction reduces the critical velocity of the vortex core 
reversal and facilitates the dynamical reversal process, 
which was demonstrated experimentally in Q . 
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